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O ■ Abstract 



We develop the minimal model of a new leading order parameterization of GPDs introduced 
by Shuvaev and Polyakov. The model for GPDs H and E is formulated in terms of the forward 
quark distributions, the Gegenbauer moments of the D-term and the forward limit of the GPD E. 
The model is designed primarely for small and medium-size values of xb, xb < 0.2. We examined 
two different models of the t-dependence of the GPDs: The factorized exponential model and 
the non-factorized Regge-motivated model. Using our model, we successfully described the DVCS 



\^ \ cross section measured by HI and ZEUS, the moments of the beam-spin ^jjj' , beam-charge A^^^'^ 



and transversely-polarized target ^jj!jf^°^''' DVCS asymmetries measured by HERMES and ^7r/ 
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measured by CLAS. The data on A^^^'^ prefers the Regge-motivated model of the t-dependence of 
D . the GPDs. The data on A^^^^'^^'^ indicates that the u and d quarks carry only a small fraction of 

X3 ■ 

the proton total angular momentum. 
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I. INTRODUCTION 



Generalized parton distributions (GPDs) parameterize non-perturbative parton correla- 
tion functions in hadronic targets l|, Q, 1^ , Q S ■ The GPDs generalize and interpolate 
between the common parton distributions and form factors. Collinear factorization theorems 
fo. deepW vi..t„a. Co^pto,, scatten.g B and fo. Ka.d elect.op.oduct.on of .nesons S give a 
theoretical possibility to experimentally constrain the GPDs. However, since the GPDs are 
functions of three arguments (excluding the known dependence on the renormalization scale) 
and since experimental observables involve convolution of the GPDs with hard scattering 
coefficients and not the GPDs themselves, it is very difficult to experimentally constrain the 
GPDs. Therefore, there is contin uing interest in modeling GPDs using various models of 



the hadronic structure 



Id, 



13, 



19|. 



The most convenient and widely used parameterization of GPDs is based on the double 



18|. Adding to the DD model 



distribution (DD) model introduced by Radyushkin 1^, 

the so-called D-teim [2^, which is required to restore the full form of polynomiality of 
the DD model, one obtains a simple, almost analytical, parameterization of GPDs, which 
can be readily used for the calculation of various observables, see e.g. 0|. However, such a 
parameterization of the GPDs has several phenomenologically unsatisfactory features. First, 
the successful description of the low Bjorken x HERA data on the total DVCS cross section 
requires a very specific (^-independent) shape for the input GPDs, which is very different 
from the input required for the description of the DVCS asymmetries measured at higher 



X 



Second, the parameterization "does not commute" with QCD evolution, i.e. it 
serves only to define the input for QCD evolution of GPDs at a certain initial scale fiQ. The 
result of the QCD evolution to the higher scale n cannot be generally expressed in the form 
used for the input. Therefore, one has to separately perform the rather non-trivial QCD 
evolution of GPDs. Third, the parameterization does not allow for an intuitive physical 



for a discussion of the physics of GPDs. 



of a new parameterization of 



motivation and interpretation, see 

In this paper, we continue and extend the analysis 24 
proton GPDs first introduced by Shuvaev and Polyakov 25[. 

Initially, the dual representation of quark GPDs of the pion was derived by Polyakov 
as a formal solution reproducing Mellin moments of the pion GPDs. In turn, the moments 
were related by crossing to the moments of the two-pion distribution amplitude, which was 
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expanded in terms of eigenf unctions of the corresponding QCD evolution equation (the 
Gegenbauer polynomials Cn"^) and eigenf unctions of the operator of the relative orbital 
moment of the pion pair (the Legendre polynomials Pi). The resulting parameterization of 
GPDs was termed dual because the idea of its derivation follows the hypothesis of duality of 
soft hadron-hadron interactions, which is the assumption that the 2 — 2 scattering arnpli- 
tude in the s-channel can be represented as an infinite series over t-channel exchanges (27 1. 
In the context of the quark GPDs of the pion, the dual parameterization implies that the 
GPDs are formally given by an infinite sum of t-channel resonances, which build up the 
two-pion distribution amplitude. 

In the successive work, Shuvaev and Polyakov postulated that a similar dual parameter- 
ization holds for proton singlet GPDs j2^ 

oo n+1 , 2\ / \ /i 

B\At.^^'')e{i-\x\)[l-'' "^''"'^ " 

n=l «=0 
odd even 

oo n+1 y 2 " 



oo R-i-i ^ 2 \ / \ 

n=l 1=0 V S / Vs / 

odd even 

oo n+1 ^ 2\ /\ /1\ 



n=l 1=0 
odd even 



In this equation, the superscript i denotes the quark flavor; x, ^, t and yU^ stand for the 



usual arguments of GPDs ^; and C^^ denote the unknown form factors. In this work, 
we restrict ourselves only to the helicity-even GPDs H and E and do not consider the two 
remaining helicity-even GPDs H and E. Since we are concerned with DVCS observables, 
which probe the singlet combinations of GPDs, we shall consider only singlet GPDs: The 
left hand side of Eq. ((H) represents the singlet combinations of the GPDs, H\x,^,t) = 
H\x, ^, t) — H\—x, ^, t) and i?*(x, ^, t) = E^x, ^, t) — E'^{—x, ^, t), which are antisymmetric 
with respect to x. 



hence, one should explain 
Eq. (P) is nothing but a 



It should be stressed that the series in Eq. diverge and_^ 
how this should be understood. As follows from the derivation 
shorthand notation for the x-moments of the GPDs. Therefore, the formal representation of 
Eq. (P) should be understood as a generalized mathematical function, which acting on the 
polynomials of x, reproduces correctly the Mellin moments of the corresponding GPDs. Note 
that the divergence in Eq. (^Q) is analogous to the divergence of any dual representation of 
the 2^2 scattering amplitude in soft hadronic physics: The s-channel series must diverge 
to reproduce infinities (poles) of the amphtude in the crossed t-channel |27|. 



Naturally, since the series in Eq. ((T)) diverge, one cannot use them to study the GPDs as 
functions of their variables. However, since Eq. (0) fixes all Mellin moments of the GPDs, 
one can readily obtain different representations of the GPDs, which would give continuous 
and finite GPDs. For instance, expandin g E g. (jT)) over the Gegenbauer polynomials Cn'^{x) 
on the interval —1 < x < 1, one obtains 



oo 



H\x, e, /i') = (1 - x') J2 t, f^')ci/\x) , (2) 



where 



n=l 
odd 



odd even 

The functions -R„p(^) are polynomials of the order n of the variable ^, which are defined in 
terms of the hypergeometric function, 

^np(0-(-l) r(^ + i)r(| + p)^ ^^^U"2'2+2 + 2'2+^'^ ) ■ 

The representation for the GPDs is obtained by replacing the form factors -B^^ by C^^. 

The dependence of the GPDs iJ* and on the renormalization scale fi is contained 
in the form factors 5^^. By construction (^, BJ^i are proportional to the rath conformal 
moment of the GPDs, which is multiplicatively renormalized to the leading accuracy j3| ■ 
Therefore, under leading order (LO) QCD evolution, the form factors have a very simple, 
well-established /x^- dependence, 

BU^^')=BM) fl^M^V"'", (5) 



Vln(/iVA2) 

where B = 11 — (2/3)A^flav (A'Wv is the number of active parton fiavors); 7„ are LO non- 
singlet anomalous dimensions [2^ 13 ■ Alternatively, as will be clear from the following 
sections, the /x^-dependence of i?^; for all / is given by the /x^-dependence of the n + 1 Mellin 
moment of the forward singlet quark distribution, J^^dx x"'q^{x, fj?). 

In addition to the simple LO evolution of 5^;, the DVCS amplitude has a very sim- 
ple form in terms of also only to the LO accuracy. Therefore, we shall use the dual 
parameterization of the GPDs as a LO parameterization. 
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II. MINIMAL MODEL OF THE DUAL PARAMETERIZATION OF GPDS W 
AND 



In this section, we consider a minimal model of the dual parameterization of the GPDs 
if* and i?*, which can be formulated in terms of the forward limit of the GPDs and 

and the Gegenbauer moments of the D-term. The t-dependence of the GPDs will be 
modelled separately. 



A. Derivation of the minimal model 



As explained in the Introduction, the dual representation in the form of Eqs. 0, © and 
(0} gives finite and continuous expressions for GPDs W and E^. However, these equations 
are impractical to use since one has to sum an infinite series of large sign-alternating terms. 
An elegant method to sum the series of Eq. (j2)) was offered by Polyakov and Shuvaev j^l ■ 
The method consists in the introduction of a set of functions whose Mellin moments generate 
the form factors and C^;, 



dxx'^Rl{x,t,i/) . 



(6) 



Note that for the singlet combinations of the GPDs that we consider in this paper, n is odd 



and k is even. Using the methods detailed in Appendix B [25|, one obtains the following 
representation of the GPDs if* in terms of the generating functions Q\ (the GPDs E'^ are 
obtained by replacing Q\ by R\) 



k=0 
even 



fc-3 



3/2 
k-l-2 



1=1 
odd 



X 



- P, - 



dyy'-'-'Ql{y,t,fj;') 



(7) 



where 

ii*W(x,e,t,/i^) 



dy 
Jo y I 



d 



X 



l-k 



1 



ds — = 
1 V2x^ 



xt 



ds— 

e{2xs-xl-e)- 



e{2xs-xl-e) 



Note that the last hne in the expression for ^'^^ is missing in the original formulation j25| : 
Its presence was noticed in Ref. j^^. 

In the present work, we consider a minimal version of the dual representation, which 
consists in retaining only the contributions of the generating functions Qq and Q2 to the 
GPDs and functions Rq and R2 to the GPDs E^. This means that we keep only the 
form factors i?^„+i, -B^„-i, C^n+i ^nn-ii we keep only the maximal I = n + 1 
and next-to-maximal I = n — 1 orbital momenta in Eq. (^Q). The main motivation for such 
an approximation is the prefactor in Eq. (jZj): In the HERA kinematics < 0.005), 
the contribution of Q\ and R\ with A; > 2 is kinematically suppressed. In the HERMES 
kinematics < 0.1), we keep Q2 and i?2 a first correction. This is also true for the 
contribution of the second line in Eq. (|7|). 

The formal representation enables one to readily establish the connection between the 
Mellin moments of the GPDs W and the form factors (similar equation holds for E^) 

odd even 

(9) 

Taking the ^ — > limit in this equation, one determines the form factors and C,^„_,_i, 



2n + 3 



dxx" (g*(x, t,/i^) + q\x,t,f/)) 



2n + 4 Jo 

CL+i{t,f^') = '^^l\xx-E\x,0,t,fi') 
2n + 3 



/ dxx"- {e\x,t, n'^) + e\x,t, jj,'^)) 
Jo 



2„ + 4, (10) 

In the forward limit t — 0, g*(x, t, ^'^) + (f'{x, t, fi^) become the singlet combination of forward 
quark parton distributions and e*(x,t,/i^) -|- (^(x, t,/i^) become the unknown forward limit 
of the singlet combination of GPDs E^. 



Since Eqs. (ITUl) fix a. 



tions Qq and Rq 



25 



m 



1 -B^„+i and C*^^^, they completely determine the generating func- 



Qlix, t, 1/) = q\x, t, I?) + q\x, t,fi^)~-J — {q'{z, t, j/) + q\z, t, /i')) , 
R},{x,t,fi^)=e'{x,t,fi^)+e'{x,t,^^)-^j^ ^ {e\z,t, fi') + e'{z,t, fi^)) . (11 



Therefore, up to the t-dependence, the functions Qq and Rq are completely constrained by 
the forward parton distributions and the forward limit of the GPDs E^. Note that Eqs. (|TT|) 
are valid at all scales fi^. 

Turning to the generating functions Q2, we notice that their modeling is more ambiguous 
as compared to the functions Qq since only the Mellin moments of Q2 are constrained. The 
constraints can be derived as follows. Considering the Mellin moments of the GPDs if*, see 
E,. O. we .o*e *e coem*,. in .o. of whicK denoted Q. i« 

.urn f V „i J, „ r (i) T{N + !)(« + !)(,, + 2) 

n=l 1=0 ^^{2'^^) ^{2 ^2) 

odd even 

On the other hand, decomposing the D-term in terms of its Gegenbauer moments 

D\z, t, /i^) = (1 - z') J2 dnit, fi') CT{z) (13) 

and using the definition 



n=l 
odd 



h'ji^l = l\zz'^D\z,t,fi'), (14) 



one obtains the desired relation between the D-term and the form factors 

n+l 

dl,{t,fi') = Y,BUt,fi')Pi{0) . (15) 

1=0 

Within the framework of the minimal version of the dual parameterization, we keep only 
the form factors with I = n±l and, hence, obtain 

This equation allows us to constrain the generating function Q\ as follows. Decomposing 
Q2 

Qlix, t, fi') = Qlix, t, i^)- — g^(^, t, + g^(a;, t, (17) 

and substituting this in Eq. |T?)|) . one obtains the following constraint on the new unknown 
function Q\ 

^'dxx-Ql{x,t,f,') = ^^^. (18) 

Turning to the generating functions R2 for the GPDs E^, we notice that since the D-term 
contribution to the GPDs E^ and W are equal and opposite in sign, see e.g. then 

(19) 



Therefore, the functions R2 can be written in the form 

i?^(x,t,/i2) = i?^(x,t,/x2)- / -Rl^{z,t,fi')-Q',{x,t,fi^). (20) 

J X ^ 

Note that this representation for R2 was first suggested by Polyakov and Shuvaev, see 
Eq. (56) of 



B. Details of the minimal model 

So far we have presented rather general consideration involved in the construction of the 
minimal model of the dual parameterization of the GPDs W and E^. In the following, we 
shall discuss such details of the parameterization as the modeling of the functions Q2 and 
e* + e* and the modeling of the t-dependence of Ql and R^. 

In the following discussion of and + e*, we assume that t = and do not explicitly 
write the t-dependence. The shape of the functions Q\ is unconstrained. We assume a 
simple form for Q\, 

Ql{x, /i') = (1 - x) (a'(/i') + b\fi^) X + c\^?) x^) , (21) 

with the coefficients a*, U and c* fixed by the constraint of Eq. (jl8p evaluated for n = 1, 3 
and 5. Note that the /i^-dependence of d\ and, hence, the /i^-dependence of a\ U and c* is 
given by Eq. (0). 

The singlet combination of the first three Gegenbauer moments of the D-term was eval- 
uated in the chiral quark soliton model at the low normalization point ~ 0.6 GeV 3^ 

dl{^^l) + d'M) = Rn [(M^ifil) + M^if^D] = i?„ , 

Ri = -8, i?3 = -2.4, i?5 = -0.8, (22) 

where is the proton momentum fraction carried by the quark and antiquark of the 
fiavor i. Note that the last equality in the first line of Eq. ()22|) comes from the fact that 
at the low normalization point /iq, in the SU(2)-symmetric chiral quark soliton model, the 
u and d quarks carry the entire proton momentum. At higher /i^ of the order a few GeV^, 



the quarks carry about half o: 
in Eq. by the factor two 



the proton momentum, which reduces the numerical values 
. In our analysis, we use 



<(/.^) = i?,M«(/xg), dM)=R,,M^{^il), ci^ = 0. (23) 



The quark momentum fractions are evaluated using the leading order CTEQ5L parton 
distributions 

Since the GPDs decouple in the forward limit, the functions e* + e* in Eq. (|iup are 
unconstrained. In this work, we follow the simple model of Ref. |4], which correctly repro- 
duces the forward limit of the first moment of the GPDs E"^ and which allows to vary the 
fraction of the nucleon spin carried by quarks 

e-(a:,/x^) = ^^5(x), (24) 

where qlg^i{x, fi"^) = q\x,fi'^) — g*(x, /i^) is the valence quark distribution. Note that unlike 
(f, e* is a symmetric function of x. The coefficients Ai and Bj^ for u and d quarks are found 
from the first and second Mellin moments of the GPDs E^ 

^^^^2,_2J^(/i2)-M^(/i2) 



Buifi^) = K-2 A„(/i2) , 5,(/i2) = kd- A,(/i2) , (25) 

where J* is the contribution to the proton total angular momentum of the quark with fiavor z; 
M2 ™^ is the proton momentum fraction carried by the valence part of the quark distribution 
function; = 1.673 and = —2.033 are quark anomalous magnetic moments. We assume 
that for the strange and charm quarks, e*(a;) = (?{x) = 0. Note also that the /x^-dependence 
ofJ.isthesa,„easforMJ.seee.g.|i. 

In summary, the x and ^-dependence of the GPDs W is specified by Eqs. dH}, (fTTjl . (fH^jl 
and (jni- The GPDs E' are specified by Eqs. (fTT |) . (j^. (j^ and (j^ . 

The t-dependence of the form factors -B^; and C^^ has to be specified extra. In this work, 
we consider two models of the t-dependence. The first model assumes that all and 
and, hence, and E^ have a factorized exponential t-dependence with the /i^-dependent 
slope 

W{x, e, t, /i^) = exp H\x, U = 0, fi') , 

E\x, e, t, /i^) = exp E\x, e, t = 0, /x^) , (26) 

where 

B{fi^) = 7.6 (1 - 0.15 ln(/xV2)) GeV^ . (27) 
9 



3- 



The normal- 



The decrease of the slope with increasing /i was taken from the model 
ization of the slope was chosen in order to reproduce the result of the exponential fit to the 
t-dependence of the differential DVCS cross section measured by the HI collaboration at 
HERA for 0.1 < \t\ < 0.8 GeV^ and at jj^ = 8 GeV^, B{fi^ = 8 GeV^) = 6.02 ± 0.35 ± 0.39 
GeV~^ jssi]. Note that a factorized model of the t-dependence with the /i^-independent 
slope commutes with the QCD evolution. While this is not so in our case ()26p. effects of the 
/i^-dependence of the slope on the QCD evolution are numerically small and, thus, 

have been neglected. 

The second model of the t-dependence is much more involved: It is non-factorized and, 
hence, the t-dependence non-trivially changes with QCD evolution. Since the dual parame- 
terization of GPDs is constructed as an infinite series of t-channel exchanges, which resembles 
the representation of hadron-hadron scattering amplitudes in Regge theory, this analogy can 
serve as a guide for the t-dependence of the GPDs. In particular, we use the following Regge 
theory motivated model for and jjl 

1 ^ ^- ' .2^ 



q'{x, t, /io) - q\x, t, fio) = ql^^{x, t, /io) = ^-^^ j QIA^^ /^o) , 

q\x, t, /ig) + q'{x, t, /io) = (^^^ [Q'ix, fxl) + q\x, fxD] , 

g{x, t, fil) = (^-^^ g{x, nl) , (28) 

where q^^x.jj?) and (f^x.jj?) are quark and antiquark forward parton distributions and 
(/(x, /i^) is the gluon forward distribution. 

The model of Eq. ()28|1 is specified at some low normalization point. In this work, = 1 
GeV^. The functions g*(t) and g*(t) at higher /x^ > /Xq are obtained by LO DGLAP evolution 
at fixed t. While the gluon function g{x, t, jj?) does not enter any LO expressions for DVCS 
observables (the hand-bag approximation), one still needs to specify g{x,t,fil) for the QCD 
evolution. 

The parameters a'-^^i, a' and a'^ can be thought of as effective slopes of the corresponding 
Regge trajectories. For the valence quarks, we use a'-^^i = 1.1(1 — x) GeV~^ [M], which 
gives a good description of the nucleon elastic form factors. Numerically similar options 
for a^^i were also considered in the literature All of them correspond to 

0.9 — 1.0 GeV^^, which is the typical slope of all known meson and baryon Regge 



val 



trajectories. 
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Drawing an analogy between the parameters ol and and the slope of the Pomeron 
trajectory a^, one would expect that a' ~ ~ a'^ = 0.25 GeV~^. However, our analysis 
of the DVCS cross section in Sect. Illll shows that larger values should be taken. In this 
work, we use 

a' = 0.9 GeV-2 , a'g = 0.5 GeV-^ . (29) 
The inconsistency between the phenomenologically large values of a' and a'g and the ones 



expected on the basis of the Regge theory was discussed in Ref . ^ | . 

Since the D-term does not have a partonic interpretation, we cannot use the model of 
Eq. to constrain the t-dependence of the Gegenbauer moments d]^, and, hence, the 
t-dependence of Q\. Instead, we imploy the results of the lattice calculations of the t- 
dependence of the first moment of GPDs, which were fitted to the dipole form [4^, and use 

<'^W=C^(« = 0)|^-i^, (30) 

where Md = 1.11 ± 0.20 GeV in the continuum limit |4a]. 

Finally, the same dipole form of the t-dependence was assumed for the 5-function contri- 
bution to the functions e* and e\ see Eq. fj^ . 



C. Predictions for the GPDs W and 

The numerical predictions for the x, C, and /i^-dependence of the GPDs if" and i?" at 
t = are presented in Figs. HJ El and El (see the captions). 

In Fig. d the results of the dual parameterization are presented as solid curves. For 
comparison, the predictions of the double distribution (DD) model with the D-term added 
are given by dashed curves. In addition, for the ^ = 0.1 case, the dot-dashed curves present 
the calculation using the dual parameterization, when the contribution of Q\ is omitted. 
Therefore, the deviation of the dot-dashed from the solid curves can serve as an estimate of 
the theoretical uncertainty associated with the modeling of the function Q\. Note that for 
^ = 0.01, the dot-dashed and solid curves are indistinguishable: Only the contribution of 
Qq is important at sufficiently low ^. 

The predictions of the DD model (dashed curves) for the singlet combination of the GPDs 



11 




FIG. 1: The singlet GPD as a function of x, and fj?. The dual parameterization results (O 
[solid curves] are compared to the predictions of the DD model (|ST|) [dashed curves]. The dot- 
dashed curves correspond to the dual parameterization, when only the function Qq is retained. 
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0.1 



FIG. 2: The singlet GPD as a function of x, ^ and /i^. The curves are obtained using the 
dual parameterization summarized in Eqs. (^0]) . (|^H) and The solid curves correspond to 
J" = 0; the dashed curves correspond to J" = 0.1; the dot-dashed curves correspond to J" = 0.2. 
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FIG. 3: The singlet GPD i?" with J" = 0. The dual parameterization calculation (solid curves) is 
confronted with the DD model calculation (dashed curves). 



W were made using the standard expressions j4 , 
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dp / da [5(x - /3 - O - <^(-a; - - a^)] h{(3, a) iS^ltll 
-1 J-i+\i3\ 2 
a; 



c?^ / da [5(x - 13 - ai) - 5{x + (3 - a^j] h{p, a) 
J-i+p 



where 



and 



h{P, a) 



3 il-\P\Y-a' 
4 (1-|/3|)3 



D^(.,/.^) = (l-.^)5^<(/.^)Cr 



(31) 



(32) 



(33) 



n=l 



Note the coefficient 1/2 in the ffist Une of Eq. (j^T|) : It is required to have the correct 



normalization of the ffist moment of -f^Do- 
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The dual parameterization predictions for the singlet GPD E'^ are summarized in Fig. |21 
We would like to emphasize that the shape of the GPDs E"^ is unknown. The model for 
the forward limit of the GPDs E^ that we used, see Eqs. (|24j) and (j25|) . was chosen rather 
arbitrary and, eventually, might turn out to be wrong. In this model, the shape of E^ is 
correlated with the fraction of the proton total orbital momentum carried by the quark, J*. 
In Fig. 121 the solid curves correspond to J" = 0; the dashed curves correspond to J" = 0.1; 
the dot-dashed curves correspond to J" = 0.2. 

For an alternative model of GPDs E^, we refer the reader to the calculations within the 
framework of the chiral quark soliton model [4^. 

In Fig. 121 the dual parameterization calculation of the singlet E"^ with J" = (solid 
curves) is compared to the DD model calculation (dashed curves). The latter was performed 
using Eq. (jHTj) after the replacement of by e* and after changing the sign in front of the 
D-teim. 



III. DVCS CROSS SECTION 



The differential cross section of deeply virtual Compton scattering (DVCS) reads 



da 1 XbU 



2 



%vcs\\ (34) 



dQ^dxBdtd(l) 32 (27r)4 ^1 + 4m%xl/Q^ ' 

where xb, and y = Q^/{xbs) {^/s is the lepton-proton invariant mass) are the usual 
Bjorken variables; is the angle between the plane formed by the leptons and the plane 

n 

formed by the final photon and the final proton 5]; Tdvcs is the full DVCS amplitude. 
The bar over the DVCS amplitude squared means that we have summed over the final 
polarization and averaged over the initial polarizations of all involved particles. 

The results of high-energy DVCS measurements at HERA are usually presented in terms 



of the DVCS cross section on the photon level 



35 



42 



2x ^ f f , J da 



aDvcs(xB, Q) = f[-)Jdtd<p ^Q,^^^ , (35) 



where 

acm(l -Z/ + Z/V2) 



(36) 



is the flux of the equivalent photons Squaring the full DVCS amplitude, averaging 

over initial polarization and summing over final polarizations, one obtains the following 
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unpolarized, t-integrated DVCS cross section on the photon level 



where 



and 



aDvcs(xB,Q^) = , " 7 ,f , = / dt\Auycs{U,Q')\\ (37) 



l-ADvcs(e, t, Q')\' = \n\'{i - e) - e in*£ + nn - l^lM ^ + ) (38) 



t, Q') = e- /' dx H\x, e, t, Q') ( i-— + ^ 



X — ^ + iO X + ^ — iO 

Throughout this paper, the skewedness parameter ^ is related to the Bjorken variable as 
^ = xb/{2 — xb) 4]. Equations ^7\i and (jHH|l for the unpolarized DVCS cross section can 
be also obtained from more general expressions derived in Ref. j^. 

It is important to note that in Eq. (jH^ we used the notation of Ref. B], which differs from 
the notation of Ref. by an overall minus sign. While this is immaterial for the DVCS 
cross section, this matters for the DVCS asymmetries. 

One appealing feature of the dual parameterization of GPDs is that the convolution 



integrals in 
Qn and 



ilq. (jHUj) can be readily taken and expressed in terms of the generating functions 
23 



+ , - 25feo • (40) 



+ + + 

The high-energy HERA data on the total DVCS cross section corresponds to very small 
^, ^ < 0.005, and to small t, t < 1 GeV^. Therefore, the contribution of the GPD E to 
the DVCS cross section is negligible. Moreover, as discussed in Sect. HH at small ^, the 
contribution of Q\ with n>2 can be safely neglected. Therefore, our predictions for ctdvcs 
within the framework of the dual parameterization are made keeping only the functions Qq, 
which, up to the t-dependence, are given by the forward quark distributions. 
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Regge 




20 40 60 80 100 

[GeV^] 

FIG. 4: The total DVCS cross section ctdvcs as a function of at fixed W . The dual pa- 
rameterization predictions are given by the solid (Regge model for the t-dependence) and dashed 
(factorized exponential t-dependence) curves. The experimental points are HI jS^i] and ZEUS |43l |. 
The statistical and systematic experimental errors are added in quadrature. 



Figures m and El present our predictions for the and W dependence of ctdvcs- The 
calculations are performed using the Regge-motivated t-dependence [solid curves] and 
the factorized exponential t-dependence fl26l) [dashed curves]. The theoretical predictions 



are compared to the HI j35[ and ZEUS |43[ data. 

Note that the ZEUS data, which were taken at 1^ = 89 GeV and = g g Qgy^ 
have been rescaled to the HI values of = 82 GeV and = 8 GeV^ using the fitted 
W and dependence of the DVCS cross section measured by ZEUS, o"dvcs oc W^^^^ and 
aovcs l/(Q2)i-54 43J. 

One can see from Fig. H] that the absolute value and the dependence of the total 
DVCS cross section is reproduced well using both the non-factorized Regge-motivated (PH|) 
and factorized exponential models of the t-dependence. However, at the highest values 
of , the exponential model of the t-dependence gives somewhat larger cdvcs because of 
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W [GeV] 

FIG. 5: The total DVCS cross section cdvcs as a function of W at fixed Q^. The captions are the 
same as in Fig. 01 



the Q^-dependent slope B (j27j) . which provides a better agreement with the highest 
ZEUS point. 

Note that the parameters a' and a'^ of the Regge-motivated model of the t-dependence, 
see Eq. fl^ . were chosen such that the theoretical calculations reproduce well the absolute 
value of (Tdvcs in Fig- IH Smaller values of a' and a'g, which would be closer to a'jp, would 
give inconsistently large values of o"dvcs- 

From Fig. one can see that the absolute value and the W dependence of ctdvcs is also 
reproduced well. The HI data 3^ somewhat prefers the results of the calculation using 
the Regge-motivated t-dependence. However, large experimental errors at large values of W 
and a slight discrepancy between the HI and ZEUS data do not allow one to draw a more 
quantitative conclusion. 

In addition to the t-integrated DVCS cross section, for the first time the HI reported the 
differential DVCS cross section jssl. The dual parameterization predictions for daBvcs/dt 
as a function of t are compared to the HI data in Fig. IHl The theoretical predictions are 
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made using Eq. (|H7j) without the integration over t 



2m2 



(41) 



dt Q^^l + Am%xl/Q^' 

As one can see from Fig. for \t\ < 0.5 GeV^ both used models of the t-dependence give 
rather similar predictions and describe well the data. However, for \t\ > 0.5 GeV^, the 
exponential model corresponds to a steeper decrease of da^YCs/dt with increasing \t\ and 
allows to describe very well the highest \t\ = 0.8 GeV^ data point. The experimental errors 
on daBYcs/dt are small enough to conclude that that the Regge-motivated model of the 
t-dependence of GPDs (|^ seems to be disfavored by the large-|i:| HI data jssj ]. 
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FIG. 6: The differential DVCS cross section as a function of \t\. The dual parameterization 
predictions (solid and dashed curves, see details in the caption to Fig. |1J) are compared to the HI 



data 



35|. 



It is important to appreciate that within the framework of the dual parameterization 
of GPDs, the HERA data on DVCS cross section were described so well without any free 
parameters: We used only the forward parton distributions for the input. To be more 
precise, the parameters a' and a'g were adjusted to reproduce only the normalization of 
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the DVCS cross section: The Q^, W and t-dependencies of the cross section were then 
predicted without any further adjustments. The calculation using the factorized exponential 
t-dependence contained no free adjustable parameters since the slope -B(/i) (jTTj) has been 



experimentally measured 



35| and, hence, it could not be varied. 



The DVCS cross section can be also described using other models of GPDs. Within 
the framework of the double distribution model, the HI and ZEUS data on ctdvcs were 
successfully described with the asjTiiptotic ansatz for h{i3, a) in Eq. (ETjK which corresponds 
to the ^-independent input, ^, /Iq) = /Xq) /iQ))/2 sjl^lj. The observation 

that the DVCS cross section at high-energies (small xb) and at high can be described 
by the GPDs, whose shape at the low input scale /ig is ^-independent, can be qualitatively 
explained as follows. Under QCD evolution, a GPD at a given small x and large //^ is 
obtained using the GPD at the low input scale fi^, which is probed ^I^^o ^ Therefore, 
the small external parameter ^ can be neglected in the input GPDs 44]. 

Other theoretical approaches, which enable one to successfully describe the HERA data 
on (Tdvcs ; include the dipole formalism j3, li^ and the formalism based on the conformal 



moments of the GPDs 



gipoi' 



IV. DVCS ASYMMETRIES 

Complete expressions for various DVCS asymmetries are well-known [sj]. In this work, 
we consider the beam-spin Am, beam-charge Ac and transversely-polarized target Aut 
asymmetries. The first two asymmetries are predominantly sensitive to the GPD H, while 
the Aut asymmetry is sensitive to both GPDs H and E. Since most of the data on these 
DVCS asymmetries have come from the HERMES collaboration at DESY [3, 4^, 1^ Isil 
we will predominantly make numerical predictions for the above asymmetries 
using the dual parameterization of GPDs H and E in the HERMES kinematics. In addition, 
predictions for the Jefferson Lab kinematics will be also presented. 
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A. Beam-spin asymmetry 

Using results of Ref. one obtains the following approximate expression for the sini 
moment of the beam-charge asymmetry 



"^0,unp 



(42) 

where e = xb^n/Q', the kinematic suppression factor and the leading harmonic of 
the Bethe-Heitler amplitude squared Cq^^^ are given in |5]; Fi and F2 are the Dirac and 
Pauli proton form factors, see e.g. 0|; H and S are defined by Eq. (pOj). Equation is 
approximate because we have neglected subleading harmonics (proportional to cf jj^p and 
•^2 unp) expansion of the Bethe-Heitler amplitude squared and the DVCS amplitude 

squared in the denominator of Eq. (j42j) . 

Note that we have introduced an additional minus sign in order to take into account 
the sign difference between our notation for Ti and the notation of Ref. ||5|]. Therefore, 
in Eq. (jl2)), the plus sign corresponds to the positively charged lepton beam. Since in 
our notation ImTi < in the bulk of the considered kinematics, A^^'^ in the HERMES 
kinematics is negative. 

n 

One should point out that we use the reference frame of Ref. |[5j| , which differs from the 
frame used by the HERMES collaboration by the direction of the 2;-axis (the Trento sign 
conventions [s^). This means that ^hermes = vt — ^,|53]. Obviously, for the sin </)-moment 
of the beam-spin asymmetry, this difference in the notations is irrelevant. 

Using the dual parameterization of GPD discussed in Sect. HH and substituting it in 
Eq. (|42j) ■ we obtain the following range of predictions at the average kinematic point of the 
HERMES measurement, {xb) = 0.11, (Q^) = 2.6 GeV^ and (t) = -0.27 GeV^ 

^Lf/"^ = —0.22 ... — 0.24 , exponential t — dependence , 

j^^J' = -0.27 ... - 0.29 , Regge t - dependence . (43) 

The smaller absolute values of A^^'^ correspond to the calculation with Ju = Jd = 0; the 
larger absolute values of A^l'^'^ correspond to the calculation with = 0.3 and Jd = (the 
variation of Jd has no noticeable affect). As expected, the small range of predictions (for a 
given model of the t-dependence) indicates small sensitivity of A^^"^ to the GPD E. 
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Our theoretical calculations compare very well to the HERMES measurement 



A%''' = -0.23 ± 0.04 ± 0.03 . 



(44) 



In addition to the avera ge H ERMES kinematics, we studied the dependence of A^^^ on 



t, xb and bin by bin 
Ju = Jd = 0. 



51 



57| . Figure [7| summarizes our predictions, which are made using 




0.05 0.1 0.15 0.2 0.25 0.3 



FIG. 7: ^^f7 as a function oi t, xb and . The dual parameterization predictions with two 



models of the t-dependence in the HERMES kinematics 
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As can be seen from Figure[7l only as a function of |t| can be helpful in distinguishing 
between the exponential and Regge models of the t-dependence, provided the experimental 
uncertainties are sufficiently small. 

We also make predictions for the beam-spin asymmetry in the CLAS kinematics. For the 
2001 average kinematic point of the CLAS kinematics ^ = 4-25 GeV, (Q^) = 1.25 GeV^, 
{xb) = 0.19 and (t) = —0.19 GeV^, our predictions compare very well to the experimental 
value, 

^uf = 0.15 . . . 0.17 , exponential t — dependence , 
^Lif = 0.18 . . . 0.20 , Regge t — dependence , 

= 0.202 ± 0.028 , CLAS [58] . (45) 

The lower values of the theoretical predictions correspond to = = 0; the larger values 
correspond to Ju = 0.3 and Jd = 0. Note the sign change in A^l'^'^ when going from the 
positron beam (HERMES) to the electron beam (CLAS). 

Recently, CLAS performed dedicated measurements of DVCS and, in particular, of the 
beam-spin asymmetry with higher energies of the lepton beam and with much wider kine- 
matic coverage in Q^, xb and t. Figure |H1 presents our predictions for the t-dependence of 
Af'^'f' aiE = 5.7 GeV, = 1.5 GeV^ and xb = 0.25. 

It is important to point out that, by construction, the minimal model of the dual param- 
eterization of GPDs is designed for small values of xb, xb < 0.2. The increase of xb from 
Xb ~ 0.2 (HERMES) to xb ~ 0.3 (current CLAS) leads to the increasing role of the gener- 
ating functions Q2 and R2 (the latter plays role at large \t\, \t\ > 0.5 GeV^), which results 
in a significant model-dependence of our results. Therefore, our predictions in Fig. |H1 should 
be taken as semi- quantitative. Strictly speaking, in order to make quantitative predictions 
for the current CLAS kinematics, one has to extend the minimal model by including higher 
generating functions Ql and Rl with k > 4. 



B. Beam-charge asymmetry 

Next we turn to the beam-charge asymmetry, Ac- Using the results of Ref. H, the 
approximate expression for Ac reads [we have neglected the same terms in the denominator 
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FIG. 8: ^^{/"^ as a function of t in the CLAS kinematics. The two curves represent calculations 



with two models of the i-dependence (see the text). 



as in Eq. 



2\ 2 C^.unp + ^,unp "^^S ( 



Ac{<P) - (^) (1 + er .r"'"^ , (46) 

V y / „ '-0,unp 



where cjunp cf are given in Ref. p]. While cjunp is smaller than cf unp; is not negli- 
gibly small. In this work, we shall concentrate on the larger contribution to Ac proportional 
to cf unp, which can be singled out by considering the cos0-moment of Ac 



Fi(t) Re n{^, t, Q2) + -^F2(t) Re£{^, t, 



„BH 
"^0,unp 



(47) 

As discussed above, 0hermes = vr — 0. Therefore, ^^'^ '^hermes _ Since in the 

considered kinematics ReTi. > 0, we obtain ^|^'"?^hermes ^ Until the end of this subsection, 
we shall imply 0hermes, but we will use for brevity. 

The range of predictions for A^'^"'^ using the dual parameterization of GPDs can be 
compared to the HERMES measurements. For the 2002 HERMES average kinematic point. 
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{xb) = 0.12, (g2) = 2.8 GeV^ and {t) = -0.27 GeV^ [50^, we obtai 



cam 

^cohtp _ g g gg ^ exponeiitial t — dependence , 
j^cosip — Q iQ Q 23 ^ Regge t — dependence , 

A™'"^ = 0.11 ±0.04 ±0.03, HERMES [50]. (48) 

The lower values of A'^^^'^ correspond to the calculation with Ju = Jd = 0; the larger values 
correspond to = 0.3 and = 0. 

For the very recent 2006 HERMES average kinematic point, {xb) = 0.10, (Q^) = 2.5 
GeV^ and (t) = -0.12 GeV^ 7^, we obtain 

Ag""^ = 0.013 . . . 0.022 , exponential t - dependence , 
A^°"<^ = 0.080 . . . 0.092 , Regge t - dependence , 

A^°^<^ = 0.063 ± 0.029 ± 0.026 , HERMES [55] . (49) 

The following two features of Eqs. (jlHj) and deserve a discussion. First, the expo- 
nential model of the t-dependence predicts the values of A™*^"^, which are much smaller than 
those calculated with the Regge model of the t-dependence. The reasons for this are the 
non-trivial t-dependence of the real part of 7i (more precisely, the non-trivial cancellation 
between two contributions to the real part of 7i) and the large values of |t| involved. 

Second, predictions with the exponential model of the t-dependence are much more sen- 
sitive to the fraction of proton spin carried by the quarks, J\ 

In addition to the average po.M.e;ece. HERMES a„a„.. Hp— 

A^^*^ as a function of t. Figure IHl presents the comparison of our theoretical predictions to 
the data. It can be seen from Fig. El that the Regge model describes the data points well 
for \t\ < 0.2 GeV^ and underestimates the asymmetry for the larger \t\ ~ 0.4 GeV^. The 
exponential model of the t-dependence dramatically fails to describe the rise of A^*"^ with 
increasing \t\. Therefore, on the basis of the comparison of our theoretical predictions to 
the t-dependence of the cos </)-moment of the beam-charge asymmetry, we conclude that the 
non-factorized Regge model of the t-dependence of GPDs is preferred over the factorized 
exponential model. 

In addition to the t-dependence of A™^'^, we make predictions for the and xb- 



dependence of A'^q^'^ in the HERMES kinematics 51, 57 1 in Fig. ^1 Our theoretical predic- 
tions are made using = Jd = 0. 
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FIG. 9: A^*^*^ as a function of t. The dual parameterization predictions with two models of the 



i-dependence are compared to the HERMES data 




FIG. 10: A^"""^ function of and xb- The dual 



of the t-dependence in the HERMES kinematics 



5l|. 



parameterization predictions with two models 



C. Transversely polarized target asymmetry 

The DVCS asymmetry with the unpolarized beam and the transversely polarized target, 
AuT, is sensitive to all four GPDs of the nucleon j^. Since we are concerned with the GPDs 
H and E, we shall consider the sin cos 0-moment of Aut, where the angle if is the angle 
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between the vector of the target polarization and the hadron scattering plane in the notation 
of Ref . . The main interest in considering this DVCS observable is that it is sensitive to 
the GPD E and, hence, to the fraction of the proton total angular momentum carried by 
quarks, J\ 

n 

According to the Trento sign convention [56[, it is recommended to use different angles, 
which are used e.g. in the HERMES analysis: 0hermes = vt -0 and 0s ,hermes - 0hermes = 

Vr + V? H- Obviously,^ A'J^'"^'^ = ^J-(<^S,hermes-0hermes)cos<^hermes_ 

Using the results p, the approximate expression for A^^''"^°^'^ reads 
^sin.cos, ^ fxA ^ _^ (2 _ 2 , + 

\ y J cq yj^p 



Note that, similarly to the above considered A^l^'^ and A™*^*^, we have introduced an addi- 
tional minus sign to compensate the sign difference between our definition of Ti and S and 
the notation of Ref. ^j. 

The theoretical predictions for the t, Q"^ and x^-dependence of A^'^''^^"'^'^ using the dual 
parameterization of GPDs and Eq. ()50|) are compared to the preliminary HERMES data 5^ 
in Fig. ^2 Note that the shown error bars correspond to the statistical and systematic 
uncertainties added in quadrature and that the systematic uncertainty does not include the 
effect of the HERMES acceptance. 

The plots for the t and Q^-dependence appear to be most informative. As can be seen 
from the upper and middle panels of Fig. [TTl our theoretical calculations reproduce the data 
fairly well, except for one point. Judging by the central experimental values, one concludes 
that the data seems to prefer the scenario with Ju = Jd = 0. 

Also, both models of the t-dependence give rather close results. Therefore, in view 
of large experimental uncertainties, it is impossible to differentiate between the Regge and 
exponential models of the t-dependence of the GPDs using the current preliminary HERMES 



data nn 4 sin </3 cos 
Liaia oil s±jjrp 



In addition, we compare our predictions for A^^^"^'^ to the HERMES measurement at 
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FIG. 11: as a function of t, xb and Q^. The dual parameterization predictions with two 

models of the t-dependence are compared to the HERMES preliminary data l:.53i- The error bars 
correspond to the statistical and systematic uncertainties added in quadrature. The systematic 
uncertainty does not include the effect of the HERMES acceptance. 



the average kinematic point, (xb) = 0.095, (Q^) = 2.5 GeV^ and (-t) = 0.12 GeV^ 

—0.14 ... — 0.10 , exponential t — dependence , 



A Sin ip cos (p 

^smApcos(f) _ _g J5 _ _ _ _ 0.10 , Regge t — dependence 



A 



cos (p 
C 



-0.149 ± 0.058 ± 0.033 , prelim. HERMES [54] . 



(51) 
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The lower absolute values of A^^^'^°^'^ correspond to the calculation with Ju = Jd = 0; 
the larger values correspond to Ju = 0.3 and J a = 0. As can be seen from Eq. (j5T| . the 
agreement between our predictions and the experimental value is very good. The central 
experimental value prefers Ju = Jd = 0. 

Predictions for A^^^^"^'^ were also made within the framework of the double distribution 



model of GPDs 



The comparison of the theretical predictions [57| to the HERMES 

f 

data [53] does not allow one to make a definite conclusion about the fraction of the proton 
total angular momentum carried by the w-quark: The data on the t-dependence somewhat 
prefers = and J„ = 0.2, while the data on xb and Q^-dependence prefer J„ = 0.2 and 
Ju = 0.4. 

In conclusion of this short discussion of the fraction of the total angular momentum of 
the proton carried by quarks, we would like to mention the lattice results of the QCDSF 
collaboration: J" = 0.37 ± 0.06 and J'^ = —0.04 ± 0.04 [4^. Our analysis presented in this 
section indicates that the value of J" should be significantly smaller. 

V. CONCLUSIONS AND DISCUSSION 

In this work, we considered the new leading order (LO) dual parameterization of GPDs 
introduced by Shuvaev and Polyakov [25]. The advantages of the dual parameterization 
include simple (forward) QCD evolution of resulting GPDs and a simple expression for the 



LO DVCS amplitude. We extended the work by Guzey and Polyakov [2J] and formulated 
the minimal model of the dual parameterization of GPDs and E^, which enables one 
to relate the GPDs to the fairly well-known quantities. In particular, apart from the t- 
dependence, the GPDs can be formulated in terms of the forward quark distributions 

and Gegenbauer moments of the D-term. The GPDs can be formulated in terms of 
the unknown forward limit of the GPDs E^ and, again, the Gegenbauer moments of the 
D-teim. The price to pay for the simplicity of our dual model is that the model is designed 
for not too large xb, xb ^ 0.2. Within the considered model, the t-dependence of GPDs has 
to be modelled separately. We considered two different models of the t-dependence: The 
factorized exponential model and the non-factorized Regge-motivated model. 

We compared predictions of our model to all available data on DVCS cross section and 
asymmetries. The Q^, W and t-dependence of the DVCS cross section at high energies 
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(small ^) measured by the HI and ZEUS collaborations were successfully described by both 
models of the t-dependence. It should be stressed that our predictions for (Jdvcs are virtually 
model- independent: Only the normalization of the cross section at one kinematic point was 
fitted by appropriately choosing the effective slope parameters in the Regge-motivated model 
of the t-dependence. 

Turning to the beam-spin DVCS asymmetry, Am, we successfully described both HER- 
MES and CLAS data on A^i^jj'^ at their respective average kinematic points. We also made 
predictions for the t, and XB-dependence of A^^"^ in the HERMES kinematics bin- by-bin 
and for the ^-dependence of A^l^'^ in the CLAS kinematics with E = 5.7 GeV. We observed 
that only the t-dcpcndence of A^^'^'^ has a chance to distinguish between the two considered 
models of the t-dcpcndcnce. 

We found that within our framework, the beam-charge asymmetry Ac is the only consid- 
ered observable which distinguishes between the Regge-motivated and exponential models 
of the ^-dependence of GPDs. While the Regge-motivated model provides a reasonable de- 
scription of ^4^^"^ in the average HERMES kinematics and of the t-dependence of A^q^'^ 
measured at HERMES, the exponential model of the t-dependence fails dramatically. 

We also compared our predictions to the HERMES measurement of the DVCS asym- 
metry measured with the unpolarized beam and the transversely-polarized target, Aut- 
We obtained a fairly good description of the preliminary HERMES data on t, Q"^ and xb- 
dependence of A^^*^^"^"^ using both models of the ^-dependence. While the experimental 
uncertainties are large, the data still seems to indicate that, within our model, Ju — Jd^ 0, 
i.e. that the u and d quarks carry only a small fraction of the proton total angular momen- 
tum. 

All comparisons to the experimental values presented in this work were done taking the 
minimal version of the dual parameterization of GPDs at its face value. We did not take into 
account such potentially important effects as next-to-leading order corrections and higher 
twist effects. Their role in the context of the dual parameterization of GPDs is a subject of 
a separate analysis. 

In conclusion, the dual parameterization of GPDs presents a new model of GPDs, which, 
with a small number of model-dependent inputs, allows for a uniform description of all 
available data on DVCS cross section and asymmetries. 
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